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Abstract 



In this paper, the second in a series of two, we complete the derivation of the lowest-order wave func- 
tion of a dimensional perturbation theory (DPT) treatment for the TV-body quantum-confined system. 
Taking advantage of the symmetry of the zeroth-order configuration, we use group theoretic techniques 
and the FG matrix method from quantum chemistry to obtain analytic results for frequencies and normal 
modes. This method directly accounts for each two-body interaction, rather than an average interaction 
so that even lowest-order results include beyond-mean-field effects. It is thus appropriate for the study 
of both weakly and strongly interacting systems and the transition between them. While previous work 
has focused on energies, lowest-order wave functions yield important information such as the nature of 
excitations and expectation values of physical observables at low orders including density profiles. Higher 
orders in DPT also require as input the zeroth-order wave functions. In the earlier paper we presented 
a program for calculating the analytic normal-mode coordinates of the large-D system and illustrated 
the procedure by deriving the two simplest normal modes. In this paper we complete this analysis by 
deriving the remaining, and more complex, normal coordinates of the system. 

1 Introduction 

Quantum confined systems involving N identical particles are common in many areas of physics, including 
chemical, condensed matter and atomic physics. In earlier papers[l, 2, 3] we have developed the method of 
dimensional perturbation theory[4] at low orders to examine the energies of quantum-confined systems such 
as atoms, quantum dots and Bose-Einstein condensates, both the ground and excited states. Dimensional 
perturbation theory (DPT) has many advantages. These include the fact that the number of particles, 



N , enters into the theory as a parameter, and so it is easy to calculate results for an arbitrary number of 
particles. Also the theory is a beyond-mean-field method, directly accounting for each interaction, rather 
than some average representation of the interactions. It is therefore appropriate for confined systems of both 
weakly interacting and strongly interacting particles. In the case of a trapped gaseous atomic Bose-Einstein 
condensate (BEC) the system is typically a weakly interacting system for which a mean-field approximation 
is valid. However, if N is increased, or Feshbach resonances are exploited to increase the effective scattering 
length, or the system is squeezed in one or more directions, it will transition into a strongly-interacting 
regime where the mean-field approach breaks down. Such systems have been created in the laboratory. [5] 
Dimensional perturbation theory is equally applicable to both regimes and may be used to study the tran- 
sition from weak to strongly interacting regimes. Even at lowest order DPT includes beyond- mean- field 
effects. 

In an earlier paper [2] we presented a detailed discussion of the dimensional continuation of the A''- 
particle Schrodinger equation, the Z) — > oo equilibrium {D^) structure, and the energy through order . 
We used the FG matrix method to derive general, analytical expressions for the many-body normal-mode 
vibrational frequencies, and we gave specific analytical results for three confined A-body quantum systems: 
the A-electron atom, A-electron quantum dot, and A-atom inhomogeneous Bose-Einstein condensate with 
a repulsive hardcore potential. 

In a subsequent paper [3] we studied the A-atom inhomogeneous Bose-Einstein condensate (BEC) with 
a repulsive hardcore potential, optimizing our low-order DPT energy by fitting to low-N diffusion Monte 
Carlo data[6] and then calculated results out to large N. Results were obtained for weak, intermediate and 
strong interatomic interactions. 

Low-order many-body dimensional perturbation methods have been applied to the A-electron atom in 
Ref. [1] where analytical expressions are obtained for the ground-state energy of neutral atoms. For Z = 1 
to 127, the energies compare well to Hartree-Fock energies with a correlation correction. 

In a subsequent paper[7] (Paper I) we have begun a detailed derivation of the zeroth-order DPT wave 
function for quantum-confined systems 

With the wave function many more properties of the system become accessible. The normal-mode 
coordinates of the zeroth-order wave function reveal the nature of the excitations of the system, and the 
lowest-order wave function provides expectation values and transition matrix elements to zeroth order in 
DPT. For macroscopic quantum-confined systems, such as the BEC, the wave function is manifested in an 
explicit way as the density profile may be viewed directly in the laboratory. This is readily calculable at 
zeroth order in DPT from the zeroth-order wave function. 

Calculating energies and wave functions to higher orders in 1/D also requires as input the lowest-order 
wave function. Indeed, as we have discussed in an earlier paper[3] which considers energies and excitation 
frequencies of the BEC, low-order DPT calculations of confined interacting particles, including strongly 
interacting systems, are accurate out to N equals a few thousand particles. To obtain yet more accurate 
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energies, and to extend the calculation beyond a few thousand particles we need to go beyond these low-order 
calculations of the energies and excitation frequencies. These calculations are facilitated by Paper I and the 
present paper since, as we have already noted above, it requires the zeroth-order DPT wave function. This 
will be pursued in later papers. 

In Paper I we set up the general theory for determining the zeroth-order DPT wave function of a quantum 
confined system. For quantum-confined systems in large dimensions, the Jacobian- weighted [8] wave function 
is harmonic, the system oscillating about a configuration termed the Lewis structure with every particle 
equidistant and equiangular from every other particle. Notwithstanding its relatively simple form, the 
large-dimension, zeroth-order wave function includes beyond-mean-field efii'ects. The Lewis structure is a 
completely symmetric configuration invariant under a point group which is isomorphic to the symmetric 
group Sn- Using group theory associated with this Sn symmetry [9] and the FG matrix method familiar 
from quantum chemistry [10], there is a remarkable reduction from N{N + l)/2 possible distinct frequencies 
to only five distinct frequencies. The Sn symmetry also greatly simplifies the determination of the normal 
coordinates. There are five sets of normal coordinates, where the normal modes of a set have the same 
frequency and transform under the same irreducible representation of ^jv • Two sets of normal coordinates 
transform under the [N] irreducible representation, two sets of normal coordinates transform under the 
[N — 1, 1] irreducible representation, and one set of normal coordinates transforms under the \N — 2, 2] 
irreducible representation. 

In Paper I we illustrated this general theory, for the zeroth-order DPT wave function of a quantum 
confined system, by determining the breathing and center-of-mass modes of the [N] species (coordinates 
transforming under the same irreducible representation are said to belong to the same species). In this 
paper, we complete this analysis by determining the modes of the more complex [N — 1, 1] and {N — 2, 2] 
species. The zeroth-order density profile of the ground-state condensate is relatively simple to calculate once 
the zeroth-order wave function is known. 

Paper I and this paper are restricted to consideration of quantum confined systems with a confining 
potential of spherical symmetry. The extension of this method to systems with cylindrical symmetry is 
relatively straightforward and will be discussed in subsequent papers. This is particularly important for the 
BEC as few current laboratory traps have spherical symmetry. 

In Sections 2 through 7 we briefly review the general theory of the zeroth-order DPT wave function 
for a quantum confined system from Paper I. In contrast to Paper I where we first discussed the symmetry 
coordinates which motivated the introduction of the primitive irreducible coordinates, in this paper we lay out 
the steps involved sequentially in the order in which they are applied to a particular problem. In Section 8 we 
begin the extension to the [AT— 1, 1] and \N— 2, 2] species by discussing the primitive irreducible coordinates 
for both the angular and radial sectors of these species. Angular and radial symmetry coordinates for the 
[A''— 1, 1] and [A''— 2, 2] species are then derived by forming appropriate linear combinations of the primitive 
irreducible coordinates. In Section 9 the frequencies and normal-mode coordinates of the system are derived 
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by transforming the G and FG matrices in the internal displacement coordinate basis to the symmetry 
coordinate basis. The problem is simplified by an extraordinary degree to two 2x2 and one 1x1 matrix 
eigenvalue equations. This reduction for [N] species has already been discussed in Paper I. Here we focus 
on the remaining, more complicated [N — 1, 1] and [A'' — 2, 2] species. Prom the solution of these reduced 
eigenvalue equations, the [N — 1, 1] and [N — 2, 2] species normal coordinates are derived in Section 9.2 . 
In Section 9.3 the normal mode motions of the particles in the original internal coordinates are derived. 
Section 10 and Section 11 summarize and conclude this paper. 



2 The £)- dimensional AT-body Schrodinger Equation 

In Paper I we considered an A^-body system of particles confined by a spherically symmetric potential and 
interacting via a common two-body potential gij . The Schrodinger equation for this system in ZP-dimensional 
Cartesian coordinates is 

N N-1 N 

H^= E 9^, ^=E^, (1) 

i=l i=l j=i-\-l 



where 



hi 



2mi ^ dxl 



and Qij = Vint 



(2) 



(3) 



are the single-particle Hamiltonian and the two-body interaction potential, respectively. The operator H 
is the iJ-dimensional Hamiltonian, and Xii, is the i/*'' Cartesian component of the i*'' particle. Konf is the 
confining potential. 



2.1 The Effective S-Wave Schrodinger Equation. 

Restricting our attention to spherically symmetric {L = 0) states of the many-body system, we transform 
from Cartesian to internal coordinates. A convenient internal coordinate system for confined systems is 



D / ^ \ 

r-j = A ^a;?^ {l<i < N) and 7^ = cos{9ij) = '^Xi^Xjv j /nrj {l<i < j < N), 



(4) 



which are the £>-dimensional scalar radii Vi of the N particles from the center of the confining potential 
and the cosines 7,j of the N{N — l)/2 angles between the radial vectors. Under this coordinate change the 
eff'ective S-wa.ve Schrodinger equation in these internal coordinates becomes 




l]k — lijlik 

)N-l N 



(5) 



where (r^ )^ = (rj)^ -|- {rjf - 2 n Tj . 
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2.2 The Jacobian- Weighted Schrodinger Equation 

Dimensional perturbation theory utilizes a similarity transformation so that the kinetic energy operator 
is transformed into a sum of two types of terms, namely, derivative terms and a repulsive centrifugal-like 
term. The latter repulsive centrifugal-like term stabilizes the system against collapse in the large- -D limit 
when attractive interparticle potentials are present. The zeroth and first orders of the dimensional (1/-D) 
expansion of the similarity-transformed Schrodinger equation are then exactly soluble for any value of N. 
In the Z) — > oo limit, the derivative terms drop out, resulting in a static problem at zeroth order, while first 
order corrections correspond to simple-harmonic normal-mode oscillations about the infinite-dimensional 
structure. 

In Paper I the weight function was chosen to be the square root of the inverse of the Jacobian, J, where [8] 

J=(rir2...rjv)^-^r(^-^-i)/2 (6) 

and r is the Gramian determinant, the determinant of the matrix whose elements are 7^- (see Appendix B 
of Paper I), so that the similarity- transformed wave function and operators (i)0 are 



J2 v]/^ and (j)0 = j2 0J 2 



(7) 



Under this Jacobian transformation, a first derivative of an internal coordinate is the conjugate momentum 
to that coordinate. The matrix elements of coordinates and their derivatives between the zeroth-order 
normal- mode functions, which are involved in the development of higher-order DPT expansions, are much 
easier to calculate since the weight function in the integrals is now unity. 

Carrying out the transformation of the Schrodinger equation of Eq. (1) via Eqs. (6) and (7) , we obtain[8] : 



(8) 



where 



and 



N 



1 

2mi dn 



+ - 



i=l 

N{N-2) + {D-N 



2mirf 



(^) 



N 



1 92 



2mi dn 2m jr 



^ EE^(t^'^ ^ 



ik 



+ 



7V(Af - 2) + (D - AT - 1)2 (l^) 



Srriirf 



(9) 



N N-1 N 

i—l i—1 J— i+1 

The latter expression for (j)T is explicitly self-adjoint since the weight function, W, for the matrix elements 
is equal to unity. The similarity-transformed Hamiltonian for the energy eigenstate is = {(^i-^T + V). 
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3 Infinite- analysis: Leading order energy 

As in Paper I, we begin the perturbation analysis by regularizing the large-dimension limit of the Schrodinger 
equation by defining dimensionally-scaled variables: 

fi = ri/K{D) , E = k{D)E and H = k{D) (^^H (11) 

where n{D) is a dimension-dependent scale factor. From Eq. (9) the kinetic energy T scales in the same way 
as so the scaled version of Eq. (8) becomes 

H^=(-j—f + %i^^ = E<^, (12) 



AD) 

where barred quantities indicate that the variables are now in scaled units. The centrifugal-like term in T 
of Eq. (9) has quadratic D dependence so the scale factor k,{D) must also be quadratic in D, otherwise the 
D oo limit of the Hamiltonian would not be finite. The precise form of k{D) depends on the particular 
system and is chosen so that the result of the scaling is as simple as possible. In previous work[2] we have 
chosen k{D) — (D—1)(D — 2N— 1)/{AZ) for the S'-wave, A^-electron atom; fl iho for the TV-electron quantum 
dot where = {D — 1){D — 2N — l)/4 and the dimensionally-scaled harmonic oscillator length and trap 
frequency respectively are Iho = \/ and lu^^ — ^'^Lu^a ! ^'^'^ D'^a^o for the BEC where ttho 



and Who = -D^Who ■ The factor of n{D) in the denominator of Eq. (12) suppresses the derivative terms as D 
increases leaving behind a centrifugal-like term in an effective potential, 



Veff(r- 7;<5 = 0) = ^ ( + ^^conf (r- 7;<5 = 0)) + ^ E ^i- (r, 7; -5 = 0) , (13) 

with 5 = 1/D, in which the particles become frozen at large D. In the D oo (5^0) limit, the excited 
states collapse onto the ground state at the minimum of Vett ■ 

We assume a totally symmetric minimum characterized by the equality of all radii and angle cosines of 
the particles when Z) — > oo, i.e. 

ri = foo {l<i<N), ^ij=^^ {l<i<j<N). (14) 

In scaled units the zeroth-order (D oo) approximation for the energy becomes 

^oo =t4ff(foo,7oo; '^ = 0). (15) 

In this leading order approximation, the centrifugal-like term that appears in Vett, even for the ground state, 
is a zero-point energy contribution satisfying the minimum uncertainty principle [11]. 
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4 Normal- mode analysis and the first-order quantum energy 
correction 



At zeroth-order, the particles are frozen in a completely symmetric, high-Z) configuration which is somewhat 
analogous to the Lewis structure in atomic physics terminology. Similarly, the first-order 1/D correction 
can be viewed as small oscillations of this structure, analogous to Langmuir oscillations. To obtain the 
1/D quantum correction to the energy for large but finite values oi D, we expand about the minimum of 
the 2? — > DO effective potential. We first define a position vector, consisting of all N{N + l)/2 internal 
coordinates: 



y 



7 



where 



and 



7 



712 




1N-2,N 
\ JN-1,N ) 

Making the following substitutions for all radii and angle cosines: 

li] =7oo +^^''^7ij> 



(16) 



(17) 



(18) 



(19) 

(20) 



where 5 = 1/i? is the expansion parameter, we define a displacement vector consisting of the internal 
displacement coordinates [primed in Eqs. (19) and (20)] 



y = 



(21) 



7 



7 



where 



f' = 



and 



r'2 



(22) 



7' 




(23) 



\ 1n-i,n J 



We find 



where 



O 



51/2=0 



7V(iV+ l)/2 



(24) 



(25) 



is the number of internal coordinates and Vo is a constant which comes from the cxphcit 5 dependence in 
the centrifugal term. The first term of the 0(((5^/^)^) term defines the elements of the Hessian matrix[12] 
F of Eq. (27) below. The derivative terms in the kinetic energy are taken into account by a similar series 
expansion, beginning with a first-order term that is bilinear in d/dy', i.e. 

^ = -^'^E E^^-^^yi + O {S'/') , (26) 

where T is the derivative portion of the kinetic energy T [see Eq. (9)]. It follows from Eqs. (24) and (26) 
that G and F, both constant matrices, are defined in the first-order S = 1/D Hamiltonian as follows: 



Hi = -\dy'^Gdy, + \y'^Fy' + Vo • 



(27) 



Thus, obtaining the first-order energy correction is reduced to a harmonic problem, which is solved by 
obtaining the normal modes of the system. 

We use the FG matrix mcthod[13] to obtain the normal-mode vibrations and, thereby, the first-order 
energy correction. A derivation of the FG matrix method may be found in Appendix A of Paper I, but the 
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main results may be stated as follows. The 6*^ normal mode coordinate may be written as: 

W]b = b^y' , (28) 
where the coefficient vector b satisfies the eigenvalue equation, 

FGb=Xbb, (29) 

with the resultant secular equation 

det(FG - XI) = 0. (30) 
The coefficient vector also satisfies the normalization condition 



b'^Gb=l. (31) 



The frequencies wl are related to A by: 



Aft = u>l (32) 
while the wave function is a product of P = N{N + l)/2 harmonic oscillator wave functions 

My') = flK,{u>l/'[q%) , (33) 

6=1 

where (^^b^^[Q']bj is a one-dimensional harmonic-oscillator wave function of frequency Qb, and rib is the 
oscillator quantum number, < < oo, which counts the number of quanta in each normal mode. 

In an earlier paper [2] we solve Eqs. (30) and (32) for the frequencies. The number of roots A of Eq. (30) 
- there are P = N{N + l)/2 roots - is potentially huge. However, due to the Sn symmetry of the problem 
discussed in Sect. 5 of Paper I, there is a stunning simplification. Equation (30) has only five distinct roots, 
A^, where runs over 0~, 0+, 1~, 1+, and 2, regardless of the number of particles in the system (see 
Refs. [2] and Sect. 7.5). Thus the energy through first-order (see Eq. (35)) can be written in terms of the 
five distinct normal- mode vibrational frequencies which are related to the roots A^ of FG by 

K=4- (34) 
The energy through first-order in 5 = 1/D is then [2] 



E = E^+5 



OO ^ 



2' 

^l={o'^, i±, ">'=o 
2 } 



(35) 



where the are the vibrational quantum numbers of the normal modes with the same frequency w^. The 
quantity ci^^n^ is the occupancy of the manifold of normal modes with vibrational quantum number n^^i and 
normal mode frequency w^ . The total occupancy of the normal modes with frequency is equal to the 
multiplicity of the root A^, i.e. 

•OO 



n„=0 
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where is the multipHcity of the f/^ root. The multiphcities of the five roots are [2] 



do+ 



1 



1, 



N 



1, 



(37) 



di- = N-1, 
d2 = N{N-3)/2. 

Note that although the equation in Ref. [1] for the energy through Langmuir order is the same as Eq. (35), 
it is expressed a Uttle differently (See Ref. [14]). 

5 The Symmetry of the Large-D, AT-body Quantum Confinement 
Problem 

5.1 Q matrices in terms of simple invEiriant submatrices 

Such a high degree of degeneracy of the frequencies in large- £>, A/'-body quantum confinement problem 
indicates a high degree of symmetry. The F, G, and FG matrices, which we generically denote by Q, are 
P X P matrices. The 5jv symmetry of the Q matrices {F, G, and FG) allows us to write these matrices in 
terms of six simple submatrices which are invariant under S^. We first define the number of 7jj coordinates 
to be 



and let In he a,n N x N identity matrix, Im an M x M identity matrix, an N x N matrix of ones and 
Jm an Mx M matrix of ones. Further, we let Rhe an N x M matrix such that Rijk = ^ij + Sik , Jnm be an 
N X M matrix of ones, and JJ^m = Jmn- These matrices are invariant under interchange of the particles, 
the 5jv symmetry, and form a closed algebra (see Appendix B of Ref. [2]). 
We can then write the Q matrices as 



where the block Qp/^/ has dimension (iVx N), block Qfi^i has dimension (A'^x M), block Q^y'^r has dimension 
(M X N), and block Q^'^' has dimension (M x M). Now, as we show in Appendix B of Ref. [2], we can 
write the Q matrices as follows: 



M = N{N -l)/2, 



(38) 




(39) 



Qf'r' = {Qa — Qb)lN + Qb^N 



(40) 



Qf'^' = (Qe — Qf)R + QfJNM 
Q^'r' = [Qc — QdlR^ + QdJjfM 



(41) 



(42) 



(43) 
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It is this structure that results in the remarkable reduction from a possible P = N{N + l)/2 distinct 
frequencies to just five distinct frequencies. 

For Q = FG, the matrix that must be diagonalized in Eq. (29), Eq. (39) becomes 

/ clIn + IJn bR+JJnm \ ,,,, 

FG=\ , (44) 

\ cR +dJMN giM + hR R + lJm I 

where we have used the following abbreviations: 

a = {FG)a-{FG)b = {Fa-F(,)Ga 

b = {FG)b = Ft,Ga 

c = {FG)c-{FG)d = {F,-Ff)Ga 

d = {FG)d = FfGa 

e = {FG\-{FG)f = {F,-Ff){Gg + {N-A)Gh) (45) 

/ = {FG)f = 2F,Gh+Ff{Gg + 2{N-:^)Gh) 

g = {FG)g-2{FG)h + {FG), = {Fg-2Fh + F,){Gg-2Gh) 

h = {FG)h-{FG),=FgGh+Fh{Gg + {N-Q)Gh)-F,{Gg + {N-h)Gh) 

I = {FG),=AFhGh + F,{Gg + 2{N-A)Gh). 

The F and G matrix elements on the right-hand sides of Eq. (45) may be derived using the graph-theoretic 
techniques discussed in Appendix B of Ref . [2] . 

We also require the G matrix for the normalization condition (Eq. (31)). It has a simpler structure than 
the FG matrix, 

. ci'In 
G=\ - . ' (46) 

g'lM + h'R^R 

d' = (G)a 

g' ^ (G),-2(G), (47) 
h' = {G)h 

{G)b = 

(G)c = 

{G)d = (48) 

(G)e = 

(G)/ = 

(G). = 0. 



where 



and 
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The quantities {G)a, {G)g and {G)h depend on the choice of k{D) (See Ref. [2]). 



6 Symmetry and Normal Coordinates 

As discussed in Paper I the FG matrix is a N(N + l)/2 x N{N + l)/2 dimensional matrix (there being 
N{N+l)/2 internal coordinates), and so Eqs. (29) and (30) could have up to N{N+l)/2 distinct frequencies. 
However, as noted above, there are only five distinct frequencies. The Sn symmetry is responsible for the 
remarkable reduction from N{N + l)/2 possible distinct frequencies to five actual distinct frequencies. As 
we shall also see, the Sn symmetry greatly simplifies the determination of the normal coordinates and hence 
the solution of the large- problem. 

6.1 Symmetrized coordinates 

The Q matrices, and in particular the FG matrix, are invariants under Sn, so they do not connect subspaces 
belonging to different irreducible representations of S'jv[10]. Thus from Eqs. (28) and (29) the normal 
coordinates must transform under irreducible representations of ^jv • Since the normal coordinates will be 
linear combinations of the elements of the internal coordinate displacement vectors f' and 7' , we first look 
at the Sn transformation properties of the internal coordinates. 

The internal coordinate displacement vectors f' and t' of Eqs. (22) and (23) are basis functions which 
transform under matrix representations of Sn, and each span the corresponding carrier spaces, however these 
representations of Sn are not irreducible representations of Sn ■ 

In Sec. 6.1 of Paper I we have shown that the reducible representation under which f' transforms is 
reducible to one 1-dimensional irreducible representation labelled by the partition [A''] (the partition denotes 
a corresponding Young diagram ( = Young pattern = Young shape) of an irreducible representation (see 
Appendix C of Paper I) and one {N — l)-dimensional irreducible representation labelled by the partition 
[A'' — 1, 1]. We also showed in Sec 6.1 of Paper I that the reducible representation under which 7' transforms 
is reducible to one 1-dimensional irreducible representation labelled by the partition [N], one {N — 1)- 
dimensional irreducible representation labelled by the partition [N — 1, 1] and one N{N — 3)/2-dimensional 
irreducible representation labelled by the partition [N — 2, 2]. Thus if da is the dimensionality of the 
irreducible representation of Sn denoted by the partition a then 



d[N] 



1 



d[N-i, 1] = 



N-1 



(49) 



d[N-2, 2] — 



N{N - 3) 
2 



and so 



d[N] + d[N-i, 1]= N 



(50) 
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giving correctly the dimension of the r' vector and that 

N{N - 1) 

C^[JV] + d[N-l, 1] + d[M-2, 2] = ^ , (51) 

giving correctly the dimension of the 7' vector. 

Since the normal modes transform under irreducible representations of Sn and are composed of linear 
combinations of the elements of the internal coordinate displacement vectors f' and 7' , there will be two 
1-dimensional irreducible representations labelled by the partition [N], two {N — l)-dimensional irreducible 
representations labelled by the partition [A'' — 1, 1] and one entirely angular A''(A'' — 3)/2-dimensional ir- 
reducible representation labelled by the partition [N — 2, 2]. Thus if we look at Eq. (37) we see that the 
0^ normal modes transform under two [N] irreducible representations, the 1^ normal modes transform 
under two [N — 1, 1] irreducible representations, while the 2 normal modes transform under the [A'' — 2, 2] 
irreducible representation. 



7 Calculating the Normal Modes. 

7.1 An Outline of the Program to Calculate the Normal Modes of the System. 

Summarizing the procedure in Paper I, we determine the normal coordinates in a three-step process: 

a) . We determine the primitive irreducible coordinates by the following procedure. First, we deter- 

mine two sets of linear combinations of the elements of coordinate vector f' which transform 
under particular non-orthogonal [N] and [N — 1, 1] irreducible representations of Sn ■ Using 
these two sets of coordinates we then determine two sets of linear combinations of the elements 
of coordinate vector 7' which transform under exactly the same non-orthogonal \N] and [A''— 1, 1] 
irreducible representations of Sn as the coordinate sets in the r' sector. Finally, another set of 
linear combinations of the elements of coordinate vector 7' which transforms under a particular 
non-orthogonal [N — 2, 2] irreducible representation of Sn is determined. In particular we set 
out to find sets of coordinates transforming irreducibly under Sn which have the simplest func- 
tional form possible subject to the requirement that they transform irreducibly under Sn- We 
call these coordinates primitive irreducible coordinates. 

b) . Appropriate linear combinations within each coordinate set from item a), above are taken so 

that the results transform under orthogonal irreducible representations of Sn ■ These are then 
the symmetry coordinates of the problem[13]. Care is taken to ensure that the transformation 
from the coordinates which transform under the non-orthogonal irreducible representations of 
Sn of item a), above, to the symmetry coordinates which transform under orthogonal irreducible 
representations of Sn preserve the identity of equivalent representations in the f' and 7' sectors. 
Furthermore, we choose one of the symmetry coordinates to be just a single primitive irreducible 
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coordinate so that it has the simplest functional form possible under the requirement that it 
transforms irreducibly under Sn. The succeeding symmetry coordinate is then chosen to be 
composed of two primitive irreducible coordinates and so on. In this way the complexity of the 
functional form of the symmetry coordinates is kept to a minimum and only builds up slowly as 
more symmetry coordinates of a given species are considered. 

c). The FG matrix is expressed in the r', 7' basis. However, if we change the basis in which the 
FG matrix is expressed to the symmetry coordinates an enormous simplification occurs. The 
N{N + l)/2 X N{N + l)/2 eigenvalue equation of Eq. (29) is reduced to one 2x2 eigenvalue 
equation for the [A^] sector, N — 1 identical 2x2 eigenvalue equations for the [A^ — 1, 1] sector 
and N{N — 3)/2 identical 1x1 eigenvalue equations for the [N — 2, 2] sector. In the case of 
the 2x2 eigenvalue equations for the [A^] and [A^ — 1, 1] sectors, the 2x2 structure allows for 
the mixing in the normal coordinates of the symmetry coordinates in the f' and 7' sectors. The 
1x1 structure of the eigenvalue equations in the [A^ — 2, 2] sector reflects the fact that there are 
no [A' — 2, 2] symmetry coordinates in the f' sector for the [A' — 2, 2] symmetry coordinates in 
the 7' sector to couple with. The [A' — 2, 2] normal modes are entirely angular. 

7.2 The Primitive Irreducible Coordinate Vectors Sr' and . 

The primitive irreducible coordinates, Sf^ , of the f' sector which transform under irreducible, though non- 
orthogonal, representations of the group Sn are given by 

Sr>=Wr>f', (52) 

where 

W,, = , (53) 

\ ^r' J 

W^^^ is a 1 X N dimensional matrix and W^^ ^' is an (AT — 1) x A" dimensional matrix. Hence we can 



write 

(54) 



Likewise the primitive irreducible coordinates, S^' , of the 7' sector which transform under irreducible, 
though non-orthogonal, representations of the group Sn are given by 

S^, = W^' 7' (55) 

and so writing 

( ^ 

(56) 
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where wl^^ is a 1 x N{N - l)/2 dimensional matrix, W^^ ^' is an {N - 1) x N{N - l)/2 dimensional 

[iV— 2 21 

matrix and WL ' is an N{N - 3)/2 x N{N - l)/2 dimensional matrix, then 



( 4"' \ 

cl^-i, 1] 



(57) 



where 



N 



(58) 



7.3 The Full Primitive Irreducible Coordinate Vector, S . 



It is useful to form a full primitive irreducible coordinate vector that groups primitive irreducible coordinates 
of the same species together as follows: 



S = P 




-glN-l, 1] 
el^-l, 1] 
-^[^-2, 2] 



v 



;;[W-2, 2] 



(59) 



where the orthogonal matrix P is given by 



and 



-^[N] 



/ 1 

10 

10 

1 

1 



1] 



1] 

IJ ji/ 

-g[N-l, 1] 



(60) 



and 5'^-'^l=4T-'^l 



Now consider applying a transformation W to Eqs. (28), (29), (30) and (31), where 



(61) 



W 















(62) 



where 



Wy'. 



(63) 
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We can now relate the primitive irreducible coordinate vector S of Eq. (59) and item a), above, to the 
symmetry coordinate vector S of item b). above. 



7.4 The Transformation, U , from Primitive Irreducible Coordinates to Sym- 
metry Coordinates. 

The symmetry coordinate vector S will be related to the primitive irreducible coordinate vector S by a 
non-orthogonal linear transformation U, i.e. 

S=US, (64) 



where 



r' 



g[N-l, 1] 

1] 
2] 



g[N-l, 1] 
g[N-2, 2] 



(65) 



and 



^[N] 



S[N-1, 1] 



g[N-l, 1] 
g[N-l, 1] 



and 2] = sl^^-'' '1 . 



(66) 



In S, symmetry coordinates of the same species are grouped together. 
Defining W through the equation 

S = Wy', 



(67) 



then Eqs. (59), (61), (62), (63), (65), (66) and (67) imply that W will be related to W by the same 
transformation, i.e. 



where 



so that 



W 









w = 


uw, 





















^ 
























u = 































1] 









I 













jr[N-2, 2] 

■7' / 




r 





\ 




/ 












T 













^[N-l, 1] 
r' 













1] 










1] 


















V 









j^^-2,21^[Ar-2,21 



(68) 



(69) 



(70) 
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Thus we have 

W^,=U^,W"^,, (71) 

where X' is f' or 7' (only 7' for the [N — 2, 2] sector). According to item b). of Sec. 7.1 W is an orthogonal 
transformation, and so 

WUW^'V = Io., (72) 

where Jq, is the unit matrix. 

Equation (72) is the essential equation for W^, to satisfy. 

From Eqs. (21), (65), (67) and (70), one derives 

SI' = W^, X' , (73) 

where X' is either of the f' or 7' internal displacement coordinate vectors and a ^[N — 2, 2] when X' = f' . 
According to Eqs. (71) and (72), U^' must satisfy 

U^, {wUWx'f} [Ux'f = la ■ (74) 

Thus to determine U^, we need to know beforehand. 

In accordance with item a), of Sec. 7.1, S^/, where a is [N] or [A'' — 1, 1], transforms under exactly the 
same representation of Sn as S^, . Likewise, according to item b). of Sec. 7.1, S^, transforms under exactly 
the same representation of Sn as S^, . Thus we must have 

U^, = U^, , (75) 

where a is [N] ov [N — 1, 1] and is a number (see below). Note that Eq. (75) substituted in Eq. (74) 
yields: 

w';,w"^,r = j^.wiwir (76) 

when a is [N\ or [N - 1, 1] . 

7.4.1 The Motion Associated with the Symmetry Coordinates about the Lewis Structure 
Configuration. 

Inverting Eq. (67) we get 

y' = W^S, (77) 

where y' is defined in Eqs. (21), (22) and (23), S is given by Eqs. (65) and (66), and W is given by Eq. (70). 
Equation (77) may be directly derived from Eq. (67) using the orthogonality of the W matrix. 
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Prom Eq. (77) we can derive 



y = 



7 



(78) 



where 



r:'a 



while 



and 



= 1 V l<i<N 



[lr]ij = i V l<i,j<N. 



(79) 
(80) 

(81) 

(82) 



Equations (78), (79), (80), (81) and (82) express the the motion associated with the symmetry coordinate 
about the Lewis structure configuration. 

7.5 The Reduction of the Eigensystem Equations in the Symmetry Coordinate 
Basis . 

7.5.1 The Central Theorem. 

Defining Qw to be a Q matrix in the symmetry coordinate basis, i.e. 

Qw = WQW'^ , 



(83) 



where Q is F, G or FG, and Qw is -FW, Gw or {FG)w ■ Since F, G and FG are invariant matrices 
(under S'Ar)[10], then 

\ 



Q 



w 












'^[N-l, 1] ® ^[N-1, 1] 








(84) 



'^[N-2, 2] ® V-2, 2] J 

where the symbol indicates the direct product, cr|^j and 't|^_j are 2x2 dimensional matrices and 
<t|^_2 2] is a single element (a number). The matrix 7[jv] is the [A/']-sector identity matrix, simply the 
number 1, I[n—i, i] is the [N — 1, l]-sector identity matrix, the {N —1) x {N — l)-dimensional unit matrix 
and J[jv— 2, 2] is the [A'' — 2, 2]-sector identity matrix, the N{N — 3)/2 x N{N — 3)/2-dimensional unit 
matrix. The matrix elements 



(85) 
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where X[ and X2 are f' or 7' (only 7' for the [N — 2, 2] sector). The Q-matrix quadrant, Qx[x^, is 
given by Eqs. (39), (40), (41), (42) and (43). Although we are not summing over the repeated index ^ in 
Eq. (85), the [cr^Jx^.x^ are independent of the W^, row label, ^ . If, when we calculate W^, , the matrix 
element [cr^]x[,x^ turns out to depend on the W^, row label, ^ , then we know that we have made a mistake 
calculating W^, . This is a strong check on the correctness of our calculations in Sect. 8. 

7.5.2 The Reduction of the Eigenvalue Equation in the Symmetry Coordinate Beisis . 

As in Paper I, we transform the basic eigenvalue equation, Eq. (29), to the symmetry coordinate basis: 

WFW^ WGW'^ Wb = Fw Gw c^^^ = h c'-''^ , (86) 

where 

Fw = WFW'^ , Gw ^ WGW'^ (87) 

c(.b)=wb. (88) 

Using Eq. (28) the 6*^ normal-mode coordinate can be written under this transformation as: 

= b^W^Wy' = [c^')] ^ ■ S, (89) 

where [q']b is now directly expressed in terms of the symmetry coordinates. 
The normal-coordinate coefficient vector, c^^' , has the form: 

/ 5„,[^]c[^l0l \ 



(90) 



where we have used Eqs. (84) and (86), and the c" satisfy the eigenvalue equations 

<7^^c« = A„c« . (91) 

The cr^>^ and cr^j^_^ arc 2 x 2-dimcnsional matrices, while <yy^_2 2] ^ one-dimensional matrix. Thus 
there are five solutions to Eq. (91) denoted by {A|^j, c^±^} , {A|^_-^ ^, '' and {A[7v-2. 2], c^^^^' ^1} , 
where X' (r' or 7' , only 7' for the [N — 2, 2] sector) labels the rows of the elements of the column vector 
c". The normal-coordinate label, b, has been replaced by the labels a, ^ and ± on the right hand side of 
Eq. (90), while the elements of the column vectors l^^ are 

[11], = Sir,, (92) 

with 1 < Tj < N -I when a=[N-l, 1], or 1 < 77 < iV(iV-3)/2 when a = [N-2, 2] . The vectors 
for the [N] and [A^ — 1, 1] sectors determine the amount of angular-radial mixing between the symmetry 
coordinates in a normal coordinate of a particular a . For the [N — 2, 2] sector, the symmetry coordinates 
are the [N — 2, 2] sector normal coordinates up to a normalization constant, [c^^~'^' ^l]:y' (see Sec. 7.5.3). 
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Writing c" as 

c« = c«xc^, (93) 
where c" is a vector satisfying the normalization condition 

\c^fc°^ = l, (94) 

and c" is a normalization factor ensuring that Eq. (31) is satisfied (see Sec. 7.5.3 below). Then the reduced 
eigenvalue equation, Eq. (91), determines c" alone. 
With the [N - 2, 2] sector, Eq. (94) yields 

c[J^^2] = 1 ^ (95) 

and so Eq. (91) yields directly 

A[iV-2, 2] = '^^^2, 2] ■ (96) 

For the [A''] and [A'^ — 1, 1] sectors, if we write 

^ / cos6'? \ 

= ± , (97) 

then from Eqs. (89), (90), (93) and (97) 

[q']b = 4 (cos ei [S^,k + sin 6% [S^,]^) . (98) 

Writing 



then Eq. (91) may be written as 



from which we derive the following analytic formula for the frequencies: 



(100) 



Equations (100) and (101) then yield 
7.5.3 The Normalization Condition. 

Using Eqs. (31), (87), (84), (85) and (90), the also satisfy the normalization condition 

[c"fcr^c" = 1 . (103) 
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Prom Eqs. (93) and (103): 



For the [A^ - 2, 2] sector, Eqs. (95) and (104) yield 

\N-2, 2] ^ 



(104) 



(105) 



[JV-2, 2] 

For the [N] and \N — 1, 1] sectors, Eqs. (97) and (104) yield for the normalization constant, 

1 



cos e% 



rG 



COS 0% 

sin 9% 



7.5.4 The Normal Coordinates. 

The normal-coordinate vector, q' , is given by 



,[JV-1, 1] 
[JV-2, 2] 




V 1 

cflcos<14^1 

cL^icoseL^i^ri 



cfsin^l^^l 
c[^lsin0[^l5l^l 



[AT-l, 1] 

cl|_ cos 




1] 1] _ 


1 [^'-l- 
f cV 


' sm t^^ 




1] 


[N-1, 1] 
Cl_ ' COS 




1] gW-i. 1] _ 


f c!_ 


' sm t/L 


1] ^[A^-l, 


1] 



g[JV-2, 2]^]Af-2, 2] 



(106) 



(107) 



7.5.5 The Motion Associated with the Normal Coordinates about the Lewis Structure Con- 
figuration. 

In terms of the normal coordinates, the unsealed internal coordinate vector is 



where 




(108) 



(109) 
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with If;/ and 1-^' given by Eqs. (81) and (82) respectively. The vectors -y'^" and y^'^ ^' are 



s(6'«) 



\ 



cos6'« 



(110) 



s(6'«) 



for a = [N] or [A^- 1, 1] , and 



/[JV-2, 21 



(111) 



(112) 



Equations (108), (109), (110), (111) and (112) express, in terms of the internal coordinates r and 7, the 
motion associated with the normal coordinates, q' , about the Lewis structure configuration . 

8 Symmetry Coordinates Belonging to the [AT — 1, 1] and 
[AT - 2, 2] Species. 

The symmetry coordinates of [N] species have been discussed in Sec. 7.6 of Paper 1. Here we derive the 
symmetry coordinates of [N — 1, 1] and [N — 2, 2] species. 



8.1 Symmetry Coordinates Belonging to the [AT — 1, 1] Species. 
8.1.1 The f' Sector. 

The Primitive Irreducible Coordinate. Consider the N{N—1) quantities — fj where 1 < i ^ j < N. 
It's clear that under Sn they transform into themselves. However, they are not all linearly independent. As 
we shall justiiy below, we identify the N — 1 linearly independent elements of S^^ ^' ^' as 



[S^!! ^' \ = f\- f^+i where 1 < i < - 1 , 



a subset of the N(N — V) r\ — r'^ . Obviously 

so we can restrict ourselves to the f- — f'- for which i<j. Furthermore we can write 



a-l 



a-1 



(113) 



(114) 



(115) 



K=0 



K=0 
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where we have used Eq. (113) in the last step. Thus we can write any of the A''(A'' — 1) f\ — f'j in terms of 
the N — 1 \S^^ ^' . That the N — 1 [S^^ ^' form a minimal, irreducible, linearly independent set 
transforming under the [N — 1, 1] representation is assured since no linear combination of the [S^^ ^' ^'], 



can be formed to give S^,, of Eq. (184) in Paper I which transforms under the [N] representation of ^jv • 
Thus the N — 1 [S^^ ^' lie in the {N — l)-dimensional subspace orthogonal to the one- dimensional 
subspace of the [A''] sector, which as we have shown in Subsection 6.1, is the [A'' — 1, 1] irreducible carrier 
space of Sn ■ We have thus proved what we set out to prove, that the N — 1 [S^^ ^' transform under 
the [N — 1, 1] irreducible representation of Sn ■ 

As we have adverted in item a), of Sec. 7.1, the \S^^ ^' of Eq. (113) have the simplest functional 
form possible. Any other choice will involve linear combinations of the [S^^ ^' of Eq. (113) and will 
likely have a more complex functional form. At the very best they will involve another {N — l)-dimensional 



subset of terms of the form — r'- . 



Furthermore, if we have two particles, we have only one primitive irreducible coordinate [s[,V = 
f i — . If we have three particles we have two primitive irreducible coordinates [Sf,' ]i = r[ — f2 and 

[2^ 1] 

[Sfii ]2 = '>'2 — r'^. Notice that adding one additional particle does not entail a change in the primitive 

[2^ 1] 

irreducible coordinate, [Sf,/ ]i , involving only the first two particles. Quite generally and by design, if we 
have A'' — 1 particles and then add one more particle to the system, according to Eq. (113) the first primitive 
irreducible coordinates [0^,, 
by this addition. Thus we have 



irreducible coordinates [S^^ ^' ^''ji , 1 < i < N — 2 , involving the first A'' — 1 particles only, are left unchanged 



-^lN-2., 1] 
[Ofi, \n-i 



(116) 



where S^^ ^' ^' is the primitive irreducible coordinate vector of the first A'' — 1 particles, and \S^^ ^' ^']jv-i 
is the primitive irreducible coordinate involving particles A'' and N — 1 . 



Prom Eq. (54) we can identify W^^ ^' as 



In matrix form this is 



[N-l, 1], 



6ij — Si+ij , where 1 < i < A^ — 1 and I < j < N . 



(117) 



1 


-1 

















1 


-1 

















1 


-1 

















1 


-1 

















1 


-1 



V ■ 



(118) 



Notice that if we have A'' — 1 particles and then add one more particle to the system, according to Eqs. (117) 
and (118), the first N — 2 rows of the matrix W^^ ^' are unchanged by the addition of an additional row. 



aside, that is, from the addition of an extra zero element at the end of the row, c.f. the discussion of the 
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paragraph prior to this one. To put it more suscinctly we have 

/ 












1] 












••• 


1 


-1 



(119) 



The Symmetry Coordinate To derive ^' and then the the symmetry coordinate S^^ ^' via 

and then Eq. (70) to arrive at wif^"^' ^\ The first step in 



[JV-l, 1] 



Eq. (73), we use Eq. (74) to calculate C/i, 
this process is to evaluate W^^, ^' \W^^, ^' ^']^ . Prom Eq. (117) we can readily show that 



N 



which in matrix form is 



{N-l, 1], 



= -s, 



i,k+l 



(120) 



2 


-1 











1— ' 


2 


-1 











-1 


2 


-1 











-1 


2 


-1 











-1 


2 



V 










\ 
































-1 







... -1 


2 


/ 



(121) 

Now Eqs. (74) and (120) (or Eq. (121)) do not uniquely define C/|f as we can transform ^' on 
the left by any orthogonal matrix, and the result will still satisfy Eq. (74). However, as we have noted above 
in item b). of Sec. 7.1., according to Eq. (64) the symmetry coordinates are composed of linear combinations 
of the primitive irreducible coordinates. Thus we choose the first symmetry coordinate to be proportional 
to just the first primitive irreducible coordinate, i.e. 



J r-^[W-l, 1]-, _/ _, 



(122) 



This involves only the first two particles in the simplest motion possible under the requirement that the 
symmetry coordinate transforms irreducibly under the [N — 1, 1] representation of Sn- li N = 2, [S^,' 
is the sole symmetry coordinate transforming under the [1, 1] representation of ^2 . For N > 2, we require 
the symmetry coordinate vector to satisfy 



^[AT-l, 1] 



g[N-h 1] 

[N-l, Ih 



(123) 



iV-l 



c.f. Eq. (116). Equations (122) and (123) mean that [S^!^ ^' is formed from the first i elements of the 
Sr' vector, which means that it is formed from the first i + 1 elements of the f' vector, i.e. it involves the 
motion of only the first i + 1 particles. 
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Thus from Eq. (70) we see that is a lower triangular matrix, and if we further require all of its 

non-zero elements in the lower triangle to be positive, these requirements, together with Eq. (74), uniquely 
specify U^~^' . Defining 



a-l 
m=l 

= 1 when a — /3 > 
= when a — /? < , 

and solving Eq. (74) for u'^~'^' , subject to the above conditions, yields 



(124) 



[N-h 111 



, ' y 5^, 



a- 



where 



Qi-j+i = ^ 5mj = 1 when j -i<l 
= when j — i> 1 . 



In matrix form Eq. (125) reads 



7-/[JV-l, 1] 



VTV2 
1 



\/2\/3 V2\/3 

1 2 3 

VsVi VsVi VsVi 



sfis/l \fi\fl ViVE ViVS 



VEVe VEVe VEVe VsVe VsVe 



\ 



(125) 



(126) 

















7j[Af-2, 1] 
r' 




















1 2 JV-3 


N-2 


JV-1 


VJV-lViV \/N-WN v/iV-lVJV 







(127) 



Thus from Eqs. (70), (117) and (125), 



I... ( ^mk 



(©i-fe+i — iSi+i,k) 



(128) 
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where 1 < i < A'' — 1 and 1 < k < N . In matrix form this is 

/ 



1, 1] _ 



1 

VTV2 


1 

VTV2 














1 


1 


2 











V2V3 


V2V3 


V2V3 


1 


1 


1 


3 








V3V4 


V3V4 


V3V4 


V3V4 


1 


1 


1 


1 


4 







ViV5 








1 


1 


1 


1 


1 


5 


VEVe 


V5\/6 


\/5%/6 


VEVe 


\/5\/6 


V5Ve 














1] 
















1 


1 1 


1 













Thus the symmetry coordinate ^' from Eqs. (73) and (128) 



IS 



[4' 



In matrix form this reads 



[JV-l, l]i _ 



— ir^_|_i ] , where 1 < i < A'' — 1 . 



^[JV-l, 1] ^ 



^(r-i+r--2rl) 

7^(n+r-^+r-^-3fD 

^(fi+f^+f^+f;-4f^) 

1 /Af-l N 



glN-2, 1] 
[b^, \n-1 



J 



(129) 



(130) 



(131) 



where S^~^' ^' is the symmetry coordinate vector of the first N — 1 particles and [St')'^"^' ^'Jat-i is the 
additional symmetry coordinate involving all A'' particles. We see from Eqs. (130) and (131) that the 
complexity of the motions described by the symmetry coordinates is kept to a minimum and builds up slowly 
as more particles become involved. Referring again to Eqs. (123) and (131), we again advert that adding 
another particle to the system does not cause widespread disruption to the symmetry coordinates. The 
symmetry coordinates, and the motions they describe, remain the same except for an additional symmetry 
coordinate which describes a motion involving all of the particles. 



Motions Associated with Symmetry Coordinates [S^ ^' ^']^ . 
(128), the motions associated with symmetry coordinates [S^^~^' ^\ 



According to Eqs. (78), (81), (79) and 
I in the unsealed internal displacement 
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coordinates r about the unsealed Lewis structure configuration 

Too = D^ahof'^lr, (132) 



are given by 



Thus we see that the motion associated with symmetry coordinate [S^~^' is an antisymmetric stretch 
motion about the Lewis structure configuration involving particles 1 and 2 . As ^ gets larger, the motion 
involves more particles, ^ + 1 particles, where the first ^ particles move one way while the (^ + 1)**^ moves 
in the other. Paradoxically though, as ^ increases the motion becomes more single-particle-like since the 
(^ -|- 1)**^ radius vector in Eq. (133) is weighted by the quantity ^ . 

8.1.2 The 7' Sector. 

The Primitive Irreducible Coordinate. Again according to item a), above, ^' should transforms 
under exactly the same non-orthogonal irreducible [N — 1, 1] representation of ^jv as S^^ ^' ^\ If fi is 
the unit vector from the center of the confining field to particle i , then fi — fi+i, where 1 < i < A'' — 1 , 
transform under exactly the same [A — 1, 1] irreducible representation of Sn as the primitive irreducible 
coordinate of the r' sector, [S^^ ^' , of Eq. (113). We also have that J^iLi'f^i is invariant under Sn, 
and so {fi — rj+i) • (X^jLi r?) = J2j'=i{'^i • '^i ~ • f^j) will transform under exactly the same [A — 1, 1] 
irreducible representation of Sn as the primitive irreducible coordinate of the f' sector, [S^^ ^' , of 
Eq. (113). Upon using 

lij =ri.fj, (134) 

this becomes Ej^i(7»j - 7(,:+i)j)- Since 7,, = 1, Y^^^iilt] - = Y^j^tlrj - Ej^,+i 7(»+i)j and 

upon using Eq. (20) we obtain — '}2j:^i+il[i+i)j) ■ Thus we identify the primitive irreducible 

coordinate of the [A — 1, 1] representation in the 7' sector as 

[S?-"' % =T.%-T. = {^r. - ^r-T,) . . (135) 

So what is W^^~^' ^' ? From Eqs. (184) of Paper 1, and Eqs. (73) and (113) 

JV JV 



fe,i=i 
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and since is symmetric Eq. (136) may be recast as 



^-^[iV-l, 1]^ _ 1 

Using Eqs. (185) of Paper I, and Eqs. (23) and (81) in Eq. (137) we have 

N l-l 



[4""' \=j:i: i^'""' ']'^ [i'^']' + Fi""'' \ [1.'].) f.. . (138) 

1=2 k=l 

Thus we identify 

= ([Wi^"'' \k + [Wl"!-'' %i [u,]k) , (139) 

where [W^^ ^' ^^]ik and [lr']i are given by Eqs. (117) (or Eq. (118)) and (81) respectively. Thus 

F^^"'' \kl = { iSik[lf']l + Su[lr,]k) - iSi+l,k[lf']l+Si+ij[lr>]k) } , (140) 

where 1 < fc < ^ < iV and 1 < z < iV - 1 . 

The Symmetry Coordinate. To derive W^~^' and then the symmetry coordinate sl^~^' via 
Eq. (73), we use Eqs. (75), (76) and (125) to calculate C/i^"^' ^' and then Eq. (71) to arrive at wlf 
The first step in this process is to evaluate ^' [W^ ^' ^']^ . Prom Eq. (139) we derive 

W?-'' Fl?"" = iN-2) Fi"-" Fi^-^' (141) 

and from Eqs. (75) and (76): 

C/^-'^l = -^[/^-.^l, (142) 

where 17^7"^' is given by Eqs. (125) (or Eq. (127)). 
Thus from Eqs. (71), (139) and (142): 

%, .1 = ([M^i^-^' [iH. + ri^-^- %i [iH.) , (143) 

where Ff'^~^' ^^'^ [^f'h are given by Eqs. (128) (or Eq. (129)) and (81) respectively. Upon using 
Eq. (128) in Eq. (143) we obtain 



1 / 

^ {Smk [If']/ +^ml [i-f'jk] - «(<5i+l,fe [If']/ + ^1+1,/ [lf']fe) 



^i{i + l){N-2) 
1 



where l<A;<Z<iVandl<i<iV-l. 



(0i_fe+i [If']; +Oj-;+i [lf']fe) - i{Si+i,k['i-f']l + 5i+i, I [If ']k)^ , (144) 
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Using Eqs. (23) and (144) in Eq. (73), we obtain the symmetry coordinate 

N N 

,1 



1 



N 



1=1 

N 



Vi(> + l)(JV-2) Vfc^S 

■ i l-l i N 

.1=2 k=l k=l l=k+l 



(145) 



^i{i + l){N-2) 



N 



lk=l l=i+2 



Motions Associated with Symmetry Coordinate ^' ^' . According to Eqs. (78), (80), (82) and 
(144), the motions associated with symmetry coordinates [S^~^' ^^]^ in the unsealed internal displacement 
coordinates 7 about the unsealed Lewis structure configuration 



looir (146) 



are given by 



1 



1 



\m=l J 



(147) 

Thus we see that the motion associated with symmetry coordinate [sl^~^' is an antisymmetric bending 
motion about the Lewis structure configuration where as the angle cosines 713, 714, 715, . . . increase, 723, 724, 
725, . . . decrease and vice versa while 712, 734, 735, 745, . . . remain unchanged. Thus as for the f' sector of the 
[A'' — 1, 1] species, [S^~^' involves the motions of particles 1 and 2 moving with opposite phase to each 
other against the remaining particles. As ^ gets larger, the motion involves more particles, ^ + 1 particles, 
where the (^ + 1)*^ particle moves with opposite phase to the first ^ particles. Paradoxically though, and 
again like the f' sector of the [N —1, 1] species, as ^ increases the motion becomes more single-particle-like 
since the angle cosines involving the (^ -|- l)*'^ particle in Eq. (147) are weighted by the quantity ^ . 

8.2 Symmetry Coordinates Belonging to the [AT — 2, 2] Species. 

There is only one sector, the 7' sector, belonging to the [A'' — 2, 2] species. 

The Primitive Irreducible Coordinate. Consider the quantities 

(fj - fj) . (f^ - ri) = jik - -fti + -fji - -ijk , iov i j k I and j > i, I > k . (148) 
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There are N\/[2^{N — 4)!] such elements [15] which transform into themselves under Sn . According to 
Eq. (20) 

7ifc - lil + Ijl - Ijk = -^(7ife - 7ii + I'ji - ijk) (149) 

so that 

iik - I'ii + iji - %k = {fi - ^j') • - r/) , where i j k I and j > i, I > k (150) 

and 

f/ = y^fi . (151) 

The A^!/[2'^(iV — 4)!] {r/ — r/) . (r^' — r/) also transform into themselves under Sn ■ Wo wish to show that 
there are N{N — 3)/2 linearly independent (r/ — r/) . {f^ — 'Pi'), a particular set of which we choose to be 
the primitive irreducible coordinates of the [N — 2, 2] species. 

There are two distinct sets of linear dependencies relating the A''!/[2^(A^ — 4)!] {f/ — r/) . (f^' — r/) , 

{fi' -f/).{r^' -fi') + {fi' -fj').{r\' -f^n') = {ri' - r/) . - fi') and (152) 
{ri' - n') • {r/ - ri) + {fi - f/) . {j^ - f/) = - f/) . {fj^ - r/) . (153) 

The number A''!/[2^(A'' — 4)!] is huge when A'' is large, but the linear dependencies of Eqs. (152) and (153) 
however, drastically reduce the number to N{N — 3)/2 linearly independent terms. 
Through Eq. (152) we can write 

{fi' - fj'} . (r>' - f/) = ^ (r^' - fi^H.') • {fk' - fi') . (154) 

K=0 

Thus we only need to consider the (A'' — 1)(A'' — 2)(A^ — 3)/2 {fi — rv+i') . {fk — fi) , where i ^ j ^ k ^ I 
and I > k . 

Equation (152) may be used iteratively again to reduce the number of (r/ — TVfi') . {f^' — r/) by writing 

l-k-l 

{fi' - CPi') • {fk - fl) = ^ {fi - n+i) ■ {fk+c' - fk^c') when I < i or k > i + 1 , (155) 

f=o 

and 

(fi' - rvTi') • {fk' - fl) = 

i-k-2 

= ^ {f - fi+i') ■ {fk+c' - fk+i+c') + {fi' - CPi') • (^l' - fi+2')+ (156) 
C=o 

l-i-3 

+ ^ {fi! - TVi') • {fi+2+c' - ''Vs+c') when k < i and I > i + 1. 
C=o 

Thus we only need to consider the {f/ — fi+i) . (rjTi' — r/) where 2<l<i— loTi + 3<l<N and 
1 < i < N — 1, and {f/ — fi+i') . (rvTi' — fi+2'), where 1 < i < N — 1 , i.e. (A'' — 1)(A'' — 3) elements to 
consider. 
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Now let's apply Eq. (153) to this reduced subset. Consider the case where 2 < I < i — 1. In this case 
Eq. (153) yields the unsurprising result 



(r/ - n+i') . {f^i - f/) = (rd' - r/) . {f/ - t^i') - {r/ - r/) . (f^i - r~^i) 

= if- - fl^Ti) . (f~^i - fi') - ifi' - fi') . (rvTi' - flTi) 
= {rPi' - fi ) . (f/ - n+i') , 



(157) 



i.e. we can write all of the {r/ — JVfi') • i'f^i ~ ^i') where 2 < I < i — 1 and 3<i<A''— lin terms of the 
(fi' - f^i) . (tCi' - fi') with i + 3 < / < iV and 1 < i < - 3 . 

Thus we reduce our N\/[2^{N — 4)!] elements of the form {fi — fj) . (r/T — f";), where i ^ j k I, 
j > i and I > k, down to the {N — 2){N — 3)/2 elements (fi' — tv+i') • {fi^i — r/) with i + 3 < I < N and 
I < i < N — 5, and the TV — 3 elements of the form {r/ — JY+i') • {fi^'i — fZyi), where 2 < i < N — 2 . Hence 
we have N(N — 3)/2 linearly elements, which is exactly the dimensionality of the [N — 2, 2] representation. 
We identify these elements as the primitive irreducible coordinates of the [N — 2, 2] representation and write 



-o[JV-2, 2] 





- f2' 


) • {fa 


- fi) 


{r2 


- rs 


) • {fl 


- fi) 


(fi' 


- fh' 


) • {fl' 


-fs') 


(^2' 


- rs 


) • {fl' 


-fs') 


{rz 


- fi 


) ■ {fl 


- fb) 


(fi' 


- fh' 


)-(f5' 


-fe') 


(^2' 


- fs 


) • (fs' 


- fe') 


(fa' 


- fi' 


)-(f5' 


- fe') 


[fi 


- fh' 


) • {fa 


- fe') 


(fi' 


- fh' 


)-{fo' 


-f?') 


{r2 


- fs 


) • {fe 


- f?') 


{rs 


- fi 


) ■ {fe 


- fr') 


{fl' 


- fs 


)-{f6' 


- f) 


{rh' 


- fe 


) ■ {fi 


- fl) 


(fi' 


- fl' 


).{f7 


-fs) 



(158) 



or 



.-g[N-2, 2] 



]ij = {fi' - f~^i') . {fJTi - fj') when i + 3<j<iVandl<i<iV-3 (159) 
= {fi - n+i) . (fps' - fj') when j = i + 2 and 2 < i < A'' - 2 . (160) 

There are a myriad of ways to choose a set of N{N — 3)/2 linearly independent, primitive irreducible 
coordinates from the 7V!/[2'^(7V — 4)!] {fi — fj) . {fk' — f/) . However, the set we have chosen has some 
significant advantages as we shall see below. 
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What about wl^ ^' ^' ? Prom Eqs. (134), (159) and (160) we can write 

N N N N 

[^^' ' ]ij = y^^i^im - ^»+l,m)(4n - Sjn)f\n' -f^' = ((^im - '^»+l,m)(4n - Sjn)l'mn > (161) 

m— 1 n— 1 m— 1 n— 1 

where k = j -1 when i + 3 < j < A'' and l<z<A''-3,orfc = j- 3 when j = i + 2 and 2<i< N -2. 
Since is a symmetric matrix, Eq. (161) can be rewritten as 

N n-1 

[^^' ' ]ij = XI X] [i^irn " ^i+l,m){Skn - Sjn) + {Sin - '^i+l,n)(^fem - Sjm)] I'mn > (162) 
n— 2 m—1 

where k = j — 1 when i + 3 < j < N and l<i<A^ — 3,orA;=j — 3 when j = i + 2 and 2 < i < N — 2 . 
Thus from Eq. (73), we note that 

N 

^^'7' = EE (163) 

[A''— 2 2] 

and Eq. (162) W^y' ' may be identified as 

[JV— 2, 2] 

[W^> ' ]ij,mn = {Sim - Si+l,rn){Skn - Sjn) + {Sin - Si+l^n){5km - Sjm) , (164) 

where 2 < n < N and l<m<n — l,k = j — 1 when i + 3 < j < N and l<i<iV — 3,orfc = j— 3 when 
j = i + 2 and 2 < i < N — 2 . Equation (164) and the above conditions on k can be written as 



(165) 



fW^iy' ' hij,mn= {Sim — Si+l,m)(^{^ — SiJ-2)Sj-l^n + SiJ-2Sj-S^n — Sjn^ 

+ {Sin — Si+i^n) ^(1 ~ Sij-2)Sj-i^m + Si,j-2 ^j-3,m " Sjm^ , 

where l<i<j-2,i + 2<j<A/',2<n<7Vandl<m<n-l. 

FjY 2 21 

The Symmetry Coordinate. We derive the symmetry coordinate Si, ' ' from Eq. (64), via Eq. (74) 
to calculate C/l^~^' and Eqs. (159) and (160) for sl^ ^' ^\ The first step in this process is to evaluate 
Wl^~^' [Wl^~^' ■ Using Eq. (164) we derive 

AT n-1 

^[N-2, 2] rjjylN-2, 2],j.i ■r^ [Ar_2, 2], rwl^-^. 2]i 

Ft r]y,iY = EE[^7' kmn[H^T k/,mn 

n=2 m=l 

N n-1 

= E/ E/ [('^im — <5i+l,m)('5fcn — <5jn) + ((^in — '5i + l,n)('5fem — <5jm)] X 

X [((5i'm — '5j' + l,ni)('^'£'n — Sj'n) + {Si'n — (5i' + l, n)(<5fe'ni — iJj'm)] 

TV TV 

= EE {Sim - Si+i^m){Skn - Sjn)y- (166) 
X [{Si'm — Si' + i^m){Sk'n — Sj'n) + {Si'n — Si'+l^n){Sk'm — Sj'm)] , 

where the m = n contribution is zero since i^i + l^k^j and i' ^ i' + 1 ^ k' ^ j' . Thus 

|Wi^ ^' ^' [W^^ ^' ^']^] j/j-/ = {Sii' - Si^i'+i - Si+iA' + Si+i^i'+i){Skk' - Skj' - Sjk' + Sjj')+ 

+{Sik' — Siji ~ Si+i^k' + Si+i^j'){Ski' — — Sji' +Sj^ii+\) 

= Ri(i+l)kj,i'{i' + l)k'j'+Ri{i+l)kj,k'j'i'{i' + l)i (167) 
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where 



Rabcd,a'b'c'd' = i^aa' - Sab' - ha' + 5bb'){5cc' - Scd' - ^dc' + ^dd') , 



(168) 



k = j-1 when i + 3<j<N and l<i<A''-3,orfc = j- 3 when j = i + 2 and 2<i<N-2. Likewise, 
k' =j'-l when i' + 3 < j' < AT and 1 < i' < TV - 3 , or fc' = i' - 3 when f = i' + 2 &nd 2 < i' < N - 2 . 

2 21 

We will find Ul, ' ' of Eqs. (71) and (74) in a two step process. First of all we will introduce a 
transformation J^~'^' which block diagonalizes Eq. (167) in the j and j' indices. Subsequently we will 
then introduce a further transformation which diagonalizes the j blocks in the indices i and i! to 

yield Eq. (74). Thus we have that 



^[Af-2, 2] ^[Ar-2, 2] 



Let's work with the transformation ^' first. If we take ^' 



to be 



2, 2]j^^_,^., _ \^-,'\^. bij^x^y (1 - 5i_i,i')i' + biv Oj_j' + i 6j/_i_i(/ - 3) , 

then after some calculation we find 



(169) 



(170) 



2] |^[^-2, 2] [j_[/V-2, 2] p{7[_^-2> 2]jT 



TV ;-2 TV /'-2 



[iV-2, 2], 



"-f77l^-2. 2] r-fJ7[W-2, 2] -r 



■7 



;=4 fc=l !'=4fe' = l 

= (2(5ji/ - - Ji+i,i')(j - 3)(j - 2)(5jj' 



[Ar-2, 2]- 



(171) 



2 21 

Thus the transformation J^, ' ' has block diagonalized Eq. (167) in the j and / indices as advertised. 

2 21 

To understand what the transformation JTI, ' ' actually achieves, let's apply it to the primitive irre- 



33 



ducible coordinate ^' . Defining 



j'=4 i'=l 
N j'-2 

j'=4 i' = l 
N j'-2 

,,/ r-o[^-2, 2], 



Af j'-2 

+ ^ ^ 6ii> Qj-j'+i Qj'-i-iif - 3) [sl^ ^' 

3'=i i' = l 

i' = l i'=i+2 

./ r-cI^-2, 2] [A,_2, 2] / •/ q^ rcI^-2. 2], 

i' = l j'=i+3 

= (r/ - fvpi') . {fi + + T^' H h f^i' + n+2' + n+s' H 1- - {j - 3)f^') 

= (r/ - r^/) . I^xi ^V' - ij 3)r/ - (r/ + r^/) j , (172) 

where 1^' is defined in Eq. (81), Qi-j is defined in Eq. (124), and Eqs (159) and (160) liave been used in the 
last step. The last term in parenthesis in the last line of Eq. (172) should be compared with Eq. (130). Not 
surprisingly, it has the same form as Eq. (130), except that it includes N — 2 particles, the i*^ and {i + 1)*^ 
particles are not included. We further note that since 

(r/-f;^i').(f/ + r;:ri') = o, (173) 

Looking now at each j block in Eq. (171), it is tridiagonal in the i index and has exactly the same form 
as Eq. (121), aside from the {j — 3)(j — 2) multiplier. From Eq. (125) we see that Z^^"^' IS 

V(j-3)(j-2) 0(i+l)(j-3)(j-2) 
i' 

Qi-i'+i Sjj' , (175) 



Vi(i + l)(j-3)(i-2) 
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where Qi-j is defined in Eq. (124). Thus from Eqs. (169), (170) and (175), 

N 1-2 



i=4 /s=l 

i-2 



= E 
= E 



k=\ v/(j-3)(i-2) 

J-2 



A: 



V'*(» + l)(j-3)(.7-2) 
1 



Vi(i + l)(i-3)(i-2) 
1 



©i-fc+i (^[lf']j" j' (1 — 5k-i,i')i' + Ski' ©i-j'+i 0j"'-fc-i(j' ~ 3)^ 

[[Ir'h i'if - 1) (1 - <5i',,'-2) ei_,v+2 + i'if - 3) ei_i.+i Qj-f + l Qf-i'-l) 

i'if - 1) (1 - Si,j,-2) ei-j'+2 + i'if - 3) Oi-i^+i e,_,v+i) . 



(176) 



The symmetry coordinate sl^ ^' ^' is thus 



AT i-2 



which from Eqs. (174) and (175) is 



kl 



(177) 



1=4, fe=l 



[4' = .... ( E^' - --^i' • I E ^' - 0- - 3)^/ 1 > (178) 



Vi(i + l)(i-3)(j-2) 



E 



where 1 < i < j — 2 and i + 2 < j < N . In matrix form Si,, ' may be written 



^[JV-2, 2] 





2, 


2] 


15 




2, 


2] 


25 




2, 


2] 


35 




2, 


2] 


16 




2, 


2] 


26 




2, 


2] 


36 




2, 


2] 


46 




2, 


2] 


17 




2, 


2] 


27 




2, 


2] 


37 


[*^' 


2, 


2] 


47 




2, 


2] 


57 



<j[JV-2, 2], 
-3=^' J 18 



(179) 
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Motions Associated with Symmetry Coordinate ^' . To calculate the motions of the unsealed 
internal displacement coordinates 7 about the unsealed Lewis structure configuration 'y^ of Eq. (146) en- 

TjY 2 2I 

gendered by the symmetry coordinates [SL, ' 'j^j, we use Eqs. (78), (80), (82) and (144). The first step 
in this process is the evaluation of [wl^~^' mn • Prom Eqs. (71), (165) and (176) we find that 

JV j'-2 
j'=4 i' = l 

+i'{3' - 3) X 



(<5i'm — <5i' + l,m) 1^(1 — ^i',j'-2)Sj'-l,n + 5i',j'-2 <5j'-3,n — ^j'nj + 

+ (^i'n — ^i' + l,n) ^(1 " -2)^]' -l,m + ^i' , j' -2 <5j'-3, m " ^j'ni^ 

—==L=== ({Qi_m+l-iSi+l,m){ej-n-{j-S)6jr.)+ (180) 

VH^+l)U-3)(j-2) V 

+ (©i-n+l - iSi+l,n){&j-m - {j " 3)(5jm)j , 



c.f. Eq. (178). Thus 

2]^ _ 2]^ _ ^ r5[N-2, 2]i r^[iV-2, 2], 

^, lij )mn — \ Uj )mn — \%] V^-^' Ju,mn 

^ [s^")""'' %i ( (e^-^+i - i^i+i,„)(e,_„ - (i - 3)j,„)+ 



Vi(«+l)(j-3)(j-2)D 

+ (©i-n+l - «<^j+l,n)(©j-m - (j - 3)6'jm)j , 



where 1 < m < n < A?", and 1 < i < i - 2 and i + 2 < j < iV . 



9 The Frequencies and Normal-Mode Coordinates of the System. 
9.1 The G and FG matrices in the Symmetry- Coordinate Basis. 

We can use the W^, matrices of Eqs. (128), (144) and (180) to calculate the reduced G and FG matrix 
elements, [cr^]x;,x^ and [cr^'^]x(,x^ , of Eq. (85) for the a = [TV - 1, 1] and {N - 2, 2] species. The 
reduced matrix elements, [fj^jjxj.x^ and [ct^Jx^.x^ for the \N\ sector are calculated in Paper I. An 
outline of this calculation is to be found in Appendix A, with the following results. 

9.1.1 The [iV - 1, 1] Species. 

The Matrix Elements [<7^_i^ i]]xj,x^ • Using Eqs. (39), (40), (41), (42), (43), (46), (47), (48), (85), 
(128) and (144) we derive 

^W-i 11 = f ^ ^' ^' ^ ^ I • (182) 
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The Matrix Elements [crfj^_^^ i]]x(,x^ • Using Eqs. (39), (40), (41), (42), (43), (44), (45), (85), (99), 
(128) and (144) we derive 



<^[iV-l, 1] 



[N-1, l]!r',r> 
FG 



'[N-x, x]h',r' = VN-2 c [<T^«_i, 1]]^-,^. = {g + {N-2)h) 
9.1.2 The [iV - 2, 2] Species. 

The Matrix Element [cr^_2, 2]] • Using Eqs. (43), (46), (47), (48), (85) and (180) we derive 



''^[JV-2, 2] — ['^[N-2, 2]Jt',t' 



= 9- 



(183) 



<^[N-2, 2] - ['^[JV-2, 2]]t',t' — 9' ■ (184) 

The Matrix Element [<r^_2^ 3]] . Using Eqs. (43), (44), (45), (85) and (180) we derive 



(185) 



9.2 The Frequencies and Normal Modes. 

Using Eq. (183) in Eq. (101) wc obtain Aj'^ j^ ^ and from Eqs. (185) and (96) wc obtain A[^_2, 2] • The 
frequencies are then determined from Eq. (34) . The Aj*^! for the [N] species are discussed in Paper I. 

Likewise the normal modes <7''|'^' and q'^_^^ for the [A^] species are also calculated in Paper I (Eqs. (203) 



[N-l, 1] 



are determined from 



and (204)). In regards to the [N — 1, 1] species, the f'-j' mixing angles, 9\j. 
Eq. (102). The normalization constant c[^~^' ^1 of the reduced normal-coordinate coefficient vector, c^''^""'^' "'^1, 
of Eqs. (90) and (93) is determined from Eqs. (105) and (106). One then determines the normal mode vector, 
q' , through Eqs. (107), (130), (145) and (178). Thus we arrive at 



' i+1 



+c. 



[N-l, 1] . oW-1, 1] 

^ ^ sm ' ' 



^i{i + l){N-2) 



' i l-l i N 

i=2 fc=l fc=l l=k+l 



(186) 



N 



.fe=l 



=i+2 



,[N-1, 1], _ [N-l, 1] 



= C_ 



cos 



JAT-l, 1] 



ir, 



, [N-l, 1] . ol^-l, 1] 

+c_ ' smO_ ' ' 



■ i l-l i N 

,1=2 k=l fe=l l=k+l 



(187) 



— « 



N 



E '^'k, i+l + E '^i+l. ' 



.fc=l 



l=i+2 



,[N-2, 2], 



J-1 



= 2] 



, - ifl^i . y Jv' - (i - 3)r/ . (188) 

v^(^+i)(j-3)(j-2) v,tt ) \h > ' \ ^ ' 
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9.3 The Motions Associated with the Normal Modes. 

Prom Eqs. (108), (109), (110), (111), (112), (128), (144) and (180) 




{[N-h 1]/ 




(189) 



where Doc is given by Eq. (109). The ^ sum for the [N] species only inchidcs one term and +r[^l, ^'yl^l, 
-r^^^ and _7[''^1 are identified in Paper I. For the [N — 1, 1] species 1 < ^ < -/V — 1, and so we have 



/ [^-1. il\ - / sin6'!^ ^' ^' ,[N-i, 1]-, s 

(+r^ = -"^°Vf(^^) ,(^[iv-i, i]),[;v-i. 1] [g+ k (e,_.+i - e%M) , (190) 

C [^-1. 1]^, _ i cos 0_ /[Jv-i, 1], 

X ((©C-i+i + ©C-j+i [Ir'lij) - [If ]ij + <5€+i,j , (191) 



^ o[N-l, 1] 



For the [A'' — 2, 2] species ^ runs over the set {A;, / : V 1 < A; < Z — 2 and k + 2 < I < N} and so 

/ lN-2, 2]n _ 1 1 r ,[N-2, 2], 

~ ^k{k + l){l-3){l-2)D c[^-2> 2] [« J« 



X ((efe_i+i-fc4+i,i)(ei-,-(/-3)^i,)+ (194) 
+(efe_,+i - kSk+i,j){ei-i -{I- 3)fe)) . 



10 Summary. 



In this paper we have completed the derivation, begun in Paper I, of the lowest-order DPT S'-wavc wave 
function of a correlated quantum confined system under spherical confinement with weak, intermediate or 
strong interparticle interactions. Remarkably, it is a largely analytic approach, yielding analytic expressions 
for energies, frequencies and wave functions (also density profiles). The achievement of analytic results is 
not due to a simplified description of the two-body interaction, but rather to the high degree of symmetry 
of the zeroth-order configuration in which every particle is equidistant and equiangular from every other 
particle in a high dimensional space. This symmetry allows us to use group theoretic methods rather than 
numerical techniques to account individually for each two-body interaction. The resulting simplification is 
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remarkable resulting in analytic formulas with N as & simple parameter as well as a stunning reduction in 
the number of normal modes. Even the zeroth-order result contains correlation effects as demonstrated by 
the interparticle angle cosine of the Lewis structure, , which differs significantly from its mean-field, i.e. 
uncorrected, value of zero. 

We have obtained the normal modes and wave functions to lowest order using the FG method[13] which 
directly relates the structure of the Schrodinger equation to normal-mode coordinates describing the funda- 
mental motions of the system (Eqs. (27) - (33)). 

However, since DPT is a beyond mean-field treatment of the interacting A''-body problem, where N 
may be large and where the N{N + l)/2 interactions can be strong, this is still a formidable problem. In 
particular there are P = N{N + l)/2 normal modes, and up to P distinct frequencies to calculate. However, 
the symmetry of the Lewis structure allows an enormous simplification in the calculation of frequencies 
and normal modes. Since the Lewis structure is invariant under interchange of particles, it satisfies an Sn 
symmetry. The full i3-dimensional Hamiltonian as well as the Hamiltonian for the zeroth-order wave function 
are invariant under particle interchange, and thus invariant under Sn ■ This symmetry results, in fact, in 
only five distinct frequencies, a remarkable reduction in the number P of possible distinct frequencies. It 
also provides an immense simplification in the calculation of the normal mode coordinates. 

We note that the Sn invariance of Eq. (27) means that the F, G and FG matrices of Eqs. (27 - 
31) are invariant under Sn ■ This implies that the eigenvectors b and so the normal mode coordinates 
transform under irreducible representations of Sn- The coordinates f- and transform reducibly under ^jv; 
however using the theory of group characters one can show that the are reducible to one one- dimensional 
[N] irreducible representation of Sn and one (N-l)-dimensional [N — 1, 1] irreducible representation of 
Sn , while the 7^^- are reducible to one one-dimensional [N] irreducible representation of Sn , one (N-1)- 
dimensional [A'' — 1, 1] irreducible representation of Sn and one N(N-3)/2-dimensional [A'' — 2, 2] irreducible 
representation of Sn ■ Since the normal coordinates are linear combinations of the internal coordinates 
and (Eqs. (21) - (23) and (28)), the normal coordinate set is comprised of two one-dimensional [A''] 
irreducible representations of Sn two (N-l)-dimensional [N — 1, 1] irreducible representations of Sn and 
one entirely angular N(N-3)/2-dimensional [N — 2, 2] irreducible representation of Sn ■ 

This group theoretic information provides an immense simplification in the calculation of the normal 
mode coordinates through the use of symmetry coordinates [13]. We have determined the normal coordinates 
and distinct frequencies in a three-step process: 

a). We define sets of primitive irreducible coordinates having the simplest functional form possible 
subject to the requirement of transforming under particular non-orthogonal irreducible represen- 
tations of Sn- For the f' sector we define two sets of linear combinations of elements of the f' 
vector which transform under non-orthogonal [N] and [A^ — 1,1] irreducible representations of 
Sn ■ We then derive two sets of linear combinations of elements of the 7' vector which transform 
under exactly these same two irreducible representations of Sn- Finally we define a set of linear 
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combinations of elements of 7' which transform under a particular non-orthogonal [N — 2, 2] 
irreducible representation of Sn- 

b) . Using linear combinations within each set of primitive irreducible coordinates, we determine 

symmetry coordinates defined to transform under orthogonal irreducible representations of Sn ■ 
Care is taken to ensure that this transformation to the symmetry coordinates preserves the 
identity of equivalent representations in the f' and 7' sectors. We choose one of the symmetry 
coordinates to be a single primitive irreducible coordinate, the simplest functional form possible 
under the requirement that it transforms irreducibly under ^jv- The next symmetry coordinate is 
chosen to be composed of two primitive irreducible coordinates and so on. Thus the complexity of 
the symmetry coordinates is minimized, building up slowly as symmetry coordinates are added. 

c) . The FG matrix, originally expressed in the f' and 7' basis, is now expressed in symmetry 

coordinates resulting in a stunning simplification. The P x P eigenvalue equation of Eq. (29) 
is reduced to one 2x2 eigenvalue equation for the [N] sector, — 1 identical 2x2 eigenvalue 
equations for the [A^ — 1, 1] sector and N(N — 3)/2 identical 1x1 eigenvalue equations for the 
[N — 2, 2] sector. For the [TV] and [N — 1, 1] sectors, the 2x2 structure allows for mixing of the 
r' and 7' symmetry coordinates in the normal coordinates (see Eq. (98)). The 1x1 striicture 
of the equations in the [N — 2,2] sector reflects the absence of f' symmetry coordinates in this 
sector, i.e. the [A'' — 2, 2] normal modes are entirely angular. 

11 Conclusions. 

The increasing interest in creating systems controlled by quantum confinement is resulting in new interest 
and new demands on the iV-body techniques of quantum physics and chemistry, originally developed for 
atoms and molecules. Mean-field treatments, such as the Hartree-Fock method in atomic physics and the 
Gross-Pitaevskii method for Bose-Einstein condensates, do not include correlation effects, and therefore 
fail for systems under tight confinement or strong interaction. These new systems, with a few hundred to 
millions of particles, present serious challenges for existing iV-body methods, most of which were developed 
with small systems in mind. 

Dimensional perturbation theory directly addresses these issues. Since the perturbation parameter is the 
dimensionality of space it is not limited to weak or strong interpartic;le interac;tions. Its formulation allows 
the symmetry of the problem to be exploited to the highest degree, enabling solutions for large N systems to 
be obtained with a minimum of numerical computation. In fact, dimensional perturbation theory for such 
systems produces analytic results that are a function of N . This means essentially that results for any N 
are obtained from a single calculation[2, 3]. 

Almost all past work using dimensional perturbation theory has focused on energies with little attention 
given to the difficult task of obtaining wave functions, even at lowest order. In Paper I we began the derivation 
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of the normal mode coordinates of the lowest-order wave function. The current paper completes this work. 
With the normal-mode coordinates of the zeroth-order wave function at hand, many other properties of the 
system beyond low-order energies become accessible. For macroscopic quantum confined systems the density 
profile is an experimentally accessible observable, and to lowest order may be calculated analytically from 
the zeroth-order wave function. The fact that it is an analytic function of A'', means that the density profile 
is obtained for any A'' in a single calculation. The normal mode coordinates directly relate to the nature of 
the motions of excited states of the system which are also experimentally accessible for macroscopic quantum 
confined systems. Expectation values may be calculated, as may transition matrix elements. The fact that 
dimensional perturbation theory is a beyond-mean-field method means that such interesting properties may 
be derived for weakly, intermediate and strongly interacting systems, the latter two appearing in experiments 
with \ow-D or large- A'' systems as well as systems with strong interparticle interactions such as is the case 
with Feshbach resonances in Bose-Einstein condensates [5]. Knowing the zeroth-order wave function is a 
precurser to calculating higher-order results. 

Almost all confining potentials for quantum confined systems in the lab currently have cylindrical sym- 
metry. These include condensates confined in a cylindrical potential as well as axially symmetric quantum 
dots, two-dimensional electronic systems in a corbino disk geometry and rotating superfluid helium systems. 
One can extend the dimensional perturbation approach to handle such systems with axial, as opposed to 
spherical symmetry. 
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A Calculation of [<T^_, ^^]x[,x',, 2]h',i'^ t^J^-i, i]]^;,^^ ^nd 

2]W,i' ' reduced G and FG matrix elements in the 
symmetry-coordinate basis. 

In this appendix we use the matrices {a — [N — 1, 1] , [N — 2, 2] , X' = r' or 7') to calculate the 

reduced G and FG matrix elements, [(T^]x[,x^ and [<t^'^]x(, , using Eq. (85): 

= (W^Os [{W^'J^r , (Al) 

where Q = G or FG and ^ is a row label. Notice that [(t^]x{,x^ on the left-hand side of Eq. (Al) should 
be independent of the row label $, on the right hand side of the equation. This is a strong check on the 
correctness of our calculations. 
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A.l The Matrix Elements i]]x[,x'^ • and [o-J^^ i]]x[,x!, • 

Equations (128) and Eq. (144) read 



[Wjf ^' \k = . / ( ^mfe - iSi+i,k] = , / (©i-fe+i - «^i+i,fc) (A2) 



and 

rlN-l, 1] 



m— 1 



Vz(z + l)(7V-2) 

where 1 < fc < ^ < iV and 1<«<A''— 1. In this case the row index £, = i ■ Thus from Eqs. (46), (Al), 
(A2) and (A3) we obtain: 

AT 



N i i 

i p=i p'=i 

~/ i 2 i i 



i(j + 1) ' 

= d' (A4) 

N N l-l 
j=l ;=2 A;=l 

= (A5) 

N k-l N 

[<rf^-^,.h',^' = EEEK^"'''']-^[G^^'-'W[(H^r-^'^'n^. 

fe=2j=l i=l 

= (A6) 
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and 



r[N-l, IhT 



N k-1 N m-l 

k=2 j=l m=2 1=1 

^1 N k-1 N m-l i 

= ^' i(i + l)(N -2) 



fc=2 j=l m=2 1=1 p=l 



({Sjl + Ski) + {Sjm + 4m)) [^{Slp' + Smp') - i{Sl,i+l + ^m,i+l)) • 



p' = l 



Using Em=2 EI^T = h Em=i E(=i -5 '5m/ in Eq. (A7), and summing over m yeilds: 

^, AT fc-1 AT i 



2i(i + 1){N -2) 

^ ^ fe=2 j=l i=l p=l 



p'=l 



m=l p'=l 



p'=i ^ 



Then using Em=i(Ep'=i -^mp' - i<5m,i+i) = and Efc=2 Ej=i = 5 Efe=i Ej=i -5 Efe=i ^jfe in Eq. 
we find: 



r N 



"■[iV-l, I]]'/',-]'' 



Ai(i + l){N-2) 



'-fe,j,!=l p=l 



E ( (E "^pJ ~ '"^i+i j) + (E "^p*^ ~ iSi+i,k 

p=l 



)) 



N i 



-2 E (E - '^i+^^i) (4-5^') ( - 2) (E -^'p' - ^-^M+i) 



= 5' + 777 



p'=i 

N i 



N , i N i 

E ( ^(E '^PJ ~ + E (E '^Pfe ~ »^i+l>fc) 



^ p=l 



-2 (E ^pJ ~ '^'5'+i.i) ) ('^i' ) ( E V ~ 



p=i 



fe=i p=i 



p' = l 



N i i 
E (E '^PJ ~ (^J') (E V ~ ^ki+l) 



j,l=l p=l 



p'=l 



- h'{N-2)^.,. 
5' + (iV - 2)h' 
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Thus we obtain for the reduced G matrix in the [N —1, 1] symmetry-coordinate basis: 

["■[AT-l, l]]*",f!' =o! [flN-l, =0 

Now letting Q = FG and using Eqs. (44), (Al), (A2) and (A3) we can derive: 

N 



(AlO) 



AT 



j,k=l p=l p'=l 
j p=l p'=l 

N i N i 



j p=l 



k=l p' = l 



a 

AT TV l-l 



(All) 



[''■[iV-l, l]]r 



j=l 1=2 fe=l 



N N l-l i 



^ ' * j=l t=2 fc=l p=l 

i 

p'=i 

1 r Af J 

E (E ^Pi ~ (^(^jk + ^jl) + /) 

i,i,fe=l p=l 

i 



2i(i + i)y/ir^ 



p'=l 



, N i i 

~ ^ ( E (E '^pj ~ ^"^'+1.0 + ( E ^'^p' ~ 



p'=i 



(A12) 



Using X^^i(X^p=i 5pj — i6i+ij) = in Eq. (A12), the term involving / is zero. Summing over j then yields: 

N / i N 



['''[iV— 1, l]]r',-f' 



1 



■ iV / I N 

E ( (E ~ »'^i+l,fc) + (E '^P' ~ 

■fe,;=i ^ p=i p=i 

X e K E ^fep' ~ + (E '^'p' ~ ) 

^ p'=i p'=i ' 

(N i i N 

E(E ~ (2^) (E ^'^p' ~ ) 

fe=l p=l p'=l ^ 



(A13) 
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The term involving 2e in Eq. (A13) gives 4ei(i + 1) (see Eq. (A4) ) leaving 



[""[JV— 1, l]]r',- 



N , i i 



k,l=l ^ p=l 



p'=i 



+ 



p=i p'— 1 
p=i p'=i 

z i 

p'=i 

AT i i 



p=i 



TV 



fe=l p=l 



p'=l 



2i(i + l)ViV^ 



2i(i + l)ViV^ 



i=i p=i p'=i 
[2A/'i(i + l)-4i(i + l)] 



= \/Arr2 e. 



[AT -2] 



An analogous calculation yields 



[<^[N-l, =\/A^5. 



(A14) 



(A15) 



Finally from Eq. (44) 



(A16) 



The first two terms in the center brackets are evaluated identically to similar terms for i]W,^' 
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yield g + [N — 2)h. The term involving ZJm is zero as shown below: 

. r N fc-1 N m-1 i 



p' = l 

N , i N 



1-1 / L J V \ 

[ E ( (E - + (E ~ ] 

/c,7 = l ^ p=l p=l ^ 



2i{i + \){N -2) 

N 

-2 E (E ~ M E E ( E "^'p' ~ 



p=i p= 

N i / N m-1 i 



p' = l 

= 0, (A17) 



where we have used I]jli(Z]p=i ^pj - = • 

Thus we obtain 



Qa-f3 = i5m,/3 = 1 when a — /3 > 



m=l 



1 / ' 



p— 1 p' — 1 

i-1 



+ (E'^P" - «^i+l,n) (E - 0' - 3)(5j>„) j , (A21) 



p=l p'=l 



^FG ^ I [""^v^i, i]]r',r'— a ['''^^l, t' — V^^V 2e (A18) 

for the reduced FG matrix in the [N —1, 1] symmetry-coordinate basis. 

A. 2 The Matrix Elements [<t^_^^ 2]]=y',=y' ^^'^ \'^[n-2, 2]]t',=y' • 
Prom Eq. (180) 

[W^'"^' %j,mn = ... , ^..."^ ^, ( (Bi_„+i - i(5i+i,„)(ej_„ - (j - 3)(5j„) 

0(i-|-l)(j-3)(j-2) V 

+ (ei-„+l - i5i+l, n){Qj-m - {j - 3)Sjm)^ , (A19) 

where 2<n<A^, l<m<n-l, l<i<j -2, 4 < j < N and (Eq. (124)) 

a-l 



= when a- P<0. (A20) 
In this case the row index ^ is the pair ij . Equations (A19) and (A20) give 
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and so from Eqs. (46), (Al) and (A21) we obtain 



Af /-I AT n-1 

= 5' + E E E E '^]^.-.^' [(<-2' 2y]_,, 



Z=2 fe=l n=2 m=l 



AT n-1 



= 5~' + '^~'E(EE[<"' E(^)'-.-[(<"'' '')"]™) • 

(A22) 



n— 2 m=l 



We can show that the term involving h' is zero by performing the double sum Y^^=2 X]fe=i • 



Af i-1 



!=2 /c=l 



1 



AT l-l 



Vi(i + l)(j-3)(j-2) fe=i V p=i 



J-1 

3^ E E ( (E ^p*^ ~ fe) ( E '^p'' ~ (■^ ~ 



p=i 



p'=i 



N , r j-1 

E ( (E ^p'' ~ ^'^'+1' ( E - (i - 



2v^i(i + l)(i-3)(j-2) 



fc,i=l ^ p=l 



p=i p'=i ^ 

JV i j-1 



fe=i p=i p'=i 
Summing over k and using the fact that i < j — 2 in the last term we obtain: 

Af l-l , r Af 



(A23) 



1 



1=2 fe=l 



1=1 p=l 



i-1 

E(E'^P'^ - i5z+l,r) (E - 0' - 3)<5ji) 

p'=l 



Vi(i + l)(j-3){j-2) 



p=l p'=l p=l 

,Ar f^i r ,f \ n .1 v^W 



Then summing over / and using J2i=iiJ2p=i ^pi - i^i+hi) = along with J2i=ii'l2p'=i ^p'I - U - = 2 

one obtains 

Af l-l 



1 



1=2 k=l 



Vz(» + l)(j-3)(j-2) 



2 (^E '^P^ ~ i^i+l, - 2 (^E ~ 



p=l 



p=l 



= 0. 



(A24) 
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Thus from Eqs. (A22) and (A24) 



['^^ — 2, 2]]t',' 



9- 



(A25) 



Similarly for the matrix element [c^^g, 2]W,1' obtain 

= [w^-'' \giM + hR'R + Um] [w\r'^ 



(A26) 



The first two terms in the center bracket are evahiated identically to the analogous terms for [cr^_2 2]]7',t' 
and from Eq. (A25) yield g. The term involving Um is zero as shown below: 

JV l-l N n-1 



1=2 fe=l 71=2 m=l 



EEEE 



i(i + l)(i-3)(i-2)^^f^^_^ 

p'=i 



X i 



p'=i 



p=i 



2»(z+l)(j-3)(j-2) 



p'=i 



j-1 

,m) V" ~ (i ~ ^)^jn) + (E ^P" ~ i^i+l,n^ (E "^P'™ ~ ~ ^^^J™) 

p'=l 



X L 



n=2 m=l p=l 



E ( (E ^p*^ ~ ^) ( E ^p'' ~ (-^ ~ 

p—i p'—i 

V - i) ( E '^p'*' ~ (■^ ~ ) 

p—i p'—i ' 

fc=l ^ p=l p'=l 

n-1 i j-1 

E E (E'^p"* ~ V" - (j - 3)(5j, 

p=l p'=l 

(E'^P" ~ (E '^P'™ ~ '^•^ ~ 



(A27) 



p=i 



p'=i 
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i{i+l){j-3){j-2) 



X L 



n,m=l 



N 



Using Efc=i(Ep=i ^pk - iSi+i,k) = yields: 

w i i-i 

fe=l p=l p'=l 

H 2 ( ^^"^ ~ ( 5^ ~ ~ ^^^^'") 

^ p=i 

p— 1 p'— 1 

i j— 1 



(A28) 



m=l ^ p=l 



p'=i 



Since i<j -2, (j^l^i Spk - iSi+i^ kj (Ep'ii Vfc " (i " 3)(5jfc) = X^^^^ 6pk - i^i+i, fc holds we obtain 



i(i + l)(i-3)(i-2) 



X L 



N i 

XI (Xl '^^'^ ~ *'^*+^' 

fe=l p=l 



m 



X! (X! '^p™ ~ 

m=l p=l 

i i 



771=1 p=l 



(A29) 



Thus the reduced FG matrix element in the symmetry-coordinate basis in the [N — 2, 2] sector is: 



fy^-PG 1 
l"[Ar-2, 2]J=7',=7' 



= 9 



(A30) 
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CO 
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